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Abstract: In a viscous Bianchi type I metric of tlie Kas- 
ner form, it is well known that it is not possible to describe 
an anisotropic physical model of the universe, which satisfies 
the second law of thermodynamics and the dominant energy 
condition (DEC) in Einstein's theory of gravity. 

We examine this problem in scalar-tensor theories of grav- 
ity. In this theory we show that it is possible to describe 
the growth of entropy, keeping the thermodynamics and the 
dominant energy condition. 

PACS number(s): 98.80.Hw, 98.80.Bp 



I. INTRODUCTION 

The standard Friedmann-Robertson- Walker (FRW) 
model has attracted considerable attention in the rela- 
tivistic cosmology literature, since the FRW model is in 
accordance with the large scale spatial homogeneity and 
isotropy of the observable universe. Perhaps, one of the 
most important characteristics of this model is, as pre- 
dicted by inflation the flatness, which agrees with the 
observed cosmic microwave background radiation. 

One can consider as source for FRW models not only 
a perfect fluid, without dissipative viscous processes, 
but also a bulk viscosity, which is compatible with the 
isotropy of the universe |^,^. Therefore, the imperfect 
fluid also can be considered as a source for FRW models. 
It is shown in ref. [Q, for a "tilting" four velocity with a 
spacelike component, that FRW cosmologies, in particu- 
lar flat models, do not necessarily represent perfect fluid 
solutions. They may correspond to exact solutions of the 
field equations of viscous fluid in which bulk and shear 
viscosity are considered. In this case, the heat conduction 
vector is not zero. 

In the early universe the kinds of matter fields are un- 
certain. The presence of anisotropy at early times is a 
very natural idea to explore, as an attempt to explain, 
among other things, the local anisotropics that we ob- 
serve today in galaxies, clusters and superclusters. Thus, 



at early time, it seems appropriate to assume a geome- 
try that is more general than merely the isotropic and 
homogeneous FRW geometry. Even though the universe, 
on a large scale, seems homogeneous and isotropic at the 
present time, there are no observational data that guar- 
antee the isotropy in an era prior to the recombination. 
In fact, it is possible to begin with an anisotropic universe 
which isotropizes during its evolution by the damping 
of this anisotropy via a mechanism of viscous dissipa- 
tion. The anisotropics described above have many possi- 
ble sources: they could be associated with cosmological 
magnetic or electric fields, long-wavelength gravitational 
waves, Yang-Mills fields, axion fields in low-energy string 
theory or topological defects such as cosmic strings or 
domain walls, among others (see ref. ||] and references 
therein) . 

We should mention that, in the Einstein theory of grav- 
ity, when an anisotropic Bianchi type-I model of the Kas- 
ner form is considered, it is not possible to describe the 
growth of entropy, due to the incompatibility between the 
second law of thermodynamics and the dominant energy 
condition [^jj^ . In this paper we would like to explore this 
problem in the situation in which a more general theory 
of gravity is considered: specifically, scalar-tensor theo- 
ries. In this respect we shall consider a sort of Jordan- 
Brans-Dicke scalar field which couples to gravity in pres- 
ence of a viscous fiuid in an anisotropic Bianchi type-I 
model of the Kasner form 
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where pi, p2 and ps are three parameters that we shall 
require to be constants. The expansion factors t^^ , t^^ 
and t^^ will be determined via Einstein's field equations. 
The space is anisotropic if at least two of the three pi 
[i = 1, 2, 3) are different. 

In the next section we give the definitions of the cor- 
responding expressions that we shall use throughout this 
paper. Also, in the same section, we show that, in stan- 
dard Einstein relativity, a viscous cosmological fluid does 
not permit the Kasner metric to be anisotropic . We 
work out the general scalar - tensor theory in Sect. Ill, and 
go on to analyze specific models in Sect. IV. Whereas in 
most cases conflicts between anisotropy /viscosity and the 
dominant energy condition still turn out to be the case, 
there are a few exceptions, as exemplified in Sect. IV. C, 
where the set of given anisotropic Kasner coefficients 
leads to a thermodynamically consistent description of 
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the viscous universe. In section V we calculate the gen- 
eration of entropy for a thermodynamically consistent 
cosmological model. Finally, in section VI we conclude 
our work. 



II. VISCOUS KASNER TYPE UNIVERSE IN 
EINSTEIN THEORY 

The Kasner universe, in Einstein's theory (with cosmo- 
logical constant A = 0), refers to a vacuum cosmological 
model given by (^), where the numbers pi, p2 and p3 
satisfy the constraints 



P1+P2+ P3 



1, pI+pI 



pI 
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An anisotropic Kasner type universe can be considered 
to be filled with an ideal (nonviscous) fluid which has 
an equation of state p = p (stiff matter - the velocity of 
sound coincides with the speed of light), where p is the 
energy density and p is the isotropic pressure. In this 
case only the constraint pi + P2 + P3 — ^ is satisfied. 

When one considers a fluid with bulk and shear viscos- 
ity coefficients, i.e. a viscous fluid, in a Kasner type uni- 
verse one can obtain mathematically self-consistent solu- 
tions. Here, the mass density and the isotropic pressure 
are proportional to the bulk and shear viscosity co- 
efficients are proportional to 1/t, and the constraints (||) 
are not satisfied any more In general the energy- 
momentum tensor for a viscous fiuid is given by 

Tal3 ^[P+{P- ^0)]UaUp - (p- £,0)gaP + "^WaP, (3) 

where Wq, p, p, ^ and r/ are the fiuid's four velocity, the 
energy density, the isotropic pressure, the bulk and shear 
(or coefficient of dynamic viscosity) viscosities, respec- 
tively. The scalar expansion and the traceless shear ten- 
sor are defined hy 6 = u"^ and 



(4) 



respectively. Here ha(3 = gap — UaU/3 is the projection 
tensor. 

From the metric (|l|) and the Einstein equations we get 
for the shear viscosity the following expression 
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{1-pi -P2 - Pa)- 
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It could be shown from thermodynamic principles that 
we should impose the conditions ^ > and 77 > in 
order to have a positive entropy generation |^]. When 
these conditions are applied in an anisotropic Kasner 
type model, the parameters entering into the metric have 
to satisfy the bound Pi + P2 + P3 ^ 1 in order to deal 
with an appropriate physical model. But, as was shown 
in refs. [^j^, it is not possible to describe the growth 
of entropy in an anisotropic model of the Kasner form 



while keeping the 2nd law of thermodynamics together 
with the DEC. 

On the other hand, we require the model to satisfy 
the DEC. This becomes specified by the range —p < 
Pj < P 0, where p is the energy density and Pj (with 
J = X, y, z) are the effective momenta related to the cor- 
responding coordinate axes. Note that both p and Pj 
scale as t^^. Thus, the DEC will give some specific re- 
lations among the Kasner parameters pi {i = 1, 2, 3). In 
the following, we introduce the symbols S = pi + P2 + Pz 
and Q = pI + p\ + p1, following ref . [|| . 

The conditions Pj < p, with j = x, y, z, yield three 
inequalities, which, by adding them, reduce to just one 
inequality given by 

2S{S-1)>0. (6) 

In a similar way, from Pj > —p we get the inequality 



S> 



(7) 
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where A = 3Q/S'2 

From expression (|7|), we see that S > 0, smce 
A> With this condition on S, we obtain from 

expression (p) that necessarily S should be greater than 
one. This conclusion yields to a negative shear viscosity, 
as can be seen from expression (|^). Therefore, from this 
result, we observe that the entropy in this sort of universe 
decreases instead of increasing. 

We should note here that the DEC, specified by the 
range for Pi (i = x,y,z), —p < Pi < p, is equiv- 
alent to the range —p < p < p, where p is the 
isotropic pressure. Effectively, for a fluid with en- 
ergy density p(t) and principal pressures Pxit), Py{t) 
and Pz(t), the energy-momentum tensor becomes 



diagU 
and (t 



P^l^ 



and using Uq, — S[ 



t: 



10|. Thus, 
, we obtain 



fr 



om eqs. 



-p + 2r)al, (i ^ x,y,z), 



i) 



(8) 



(9) 



where there is no sum over i, and ^ = 0. After adding 
the inequalities Pi < p {i = x,y,z), one obtains p < p. 
On the other hand, after adding —p < Px — x, y, z) 
we get —p < p. In both cases we have used the result 
al + a2 + (J^ = 0. This result is easy generalized to the 



case where the bulk viscosity is taken into account, 
we obtain ~p + < p < p + 6£^. 



Here, 



III. SCALAR-TENSOR THEORIES AND THE 
FIELD EQUATIONS 

We start with the Lagrangian density for the scalar 
tensor theory of gravity 
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L 



cf>iR - 2A{cb)) + ^g'^-d^^ 



(10) 



where g — det{gaf3). The arbitrary functions oj — lj (0) 
and A — A (0) distinguish the different scalar-tensor 
gravitational theories: is the coupling function and 

A(0) is a potential function and plays the role of a cos- 
mological constant . 

The Euler-Lagrange equations of motion for gap and 
(f) obtained from the action 5* = J(L + Lj^,) d^x/lQ-K are 



■ 9apR ■ 



and 



Tap - K4>)9al3 



(11) 



2a; (0) + 3 



2cj(</>) 
duj{4>) 
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respectively, where the partial derivatives are denoted 
by a prime and covariant derivatives are denoted by a 
semicolon. Tq,^ is the stress energy-tensor which be- 
comes calculated from L^^ through the definition Tap = 
-g ij^] . General relativity is recovered in the 



limit OJ — > oo and <j) = const. = 1/G. We can conve- 
niently write (^l]) as 



■Q/3 



gap 
1 
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We now apply this theory to a homogeneous and 
anisotropic cosmology of the Kasner type, expressed by 
the metric (|^), together with the energy-momentum ten- 
sor specified by eq. (^). The signature here is (H ), 

and the Ricci tensor is given by Rap = Rapj- scalar 
expansion and the traceless shear tensor are 

e = sr\ 



3 -P. 



with 
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2t2 



c2- 



respectively. Then, eq. ( |T^ reduces to the following field 
equations: 



5 • 202A'(0)-20A(0) 



Stt 



2w(0) + 3 



p-3 



2w((/i) + 3 



2uj{(t))+3 



02. (14) 



Eq. (|l3| ) gives rise to the following set of eqs.: 
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and 



4-77 
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1 — b — rjt — —t 



^(0) 



= -A(0) - 



(15) 



P 



S\ 45 ■ 
7^ +3^^ 



1 

20 
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(16) 



where the dots specify derivatives with respect to the 
cosmological time. 

Since the pi (i — 1,2,3) are all linearly independent 
we get from cq. (16) that 



l-S--t] > 0. 



(17) 



From the same eq., we get a condition that, after sub- 
stituting eq. (O) into this expression we get 



20 



1 

~ 20 
3^2 



550 



A(0). 



(18) 



This latter equation, together with eq. (y_5|), allows us 
to obtain explicit expressions for p and p which result to 
be given by 



8ttp 
and 
Snp - 



(52 - Q) 50 c^(0)02 



2^2 



t 



iS-Q)<f> 
2^2 



20 

205 5(1 - 5) 



0A(0) (19) 
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6*2 
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2^ +0A(0)-|-87r-e. 



(20) 



These two expressions can be substituted into the 
equation of motion for the scalar field 0, eq. (|lj), and 
obtain 



5 
T 



202 A (0) - 20A(0) 



1 



2w(0) 



(52 + Q - 25) : 
i2 



35 
t 



(2u;(0) + 3) 
. 2 

UJ 



-4A0 



(21) 



i2 



2tj(0) -1-3 

We should note that, when = const. = 1/G and A(0) = 
0, this set of equations coincides with that considered by 
the authors of ref In the following we shall consider 
some particular cases. 



3 



IV. SOME SPECIFIC MODELS 



B. Brans-Dicke in a Kasner-type universe 



The set of equations described above can be solved for 
the particular situation to ((f)) = ujq and A(0) = Ao0~" , 
where wq, Aq and n are arbitrary real constants. This 
choice allows us to write the following particular solution 

for the scalar field: (p — cfo [ — ) as can see from (21). 

Here 0o and to are two constants that would be associated 
with the present values of </) and t, respectively. Replacing 
these expressions into eq. ( ^l|) we obtain the following 
expression for the Aq parameter: 



In this subsection we consider the Brans-Dicke theory, 
i.e. we take Aq = with n and luq for the moment ar- 
bitrary parameters. Thus, from (p3), wc have that the 
shear viscosity is still given by (|25|) and the energy den- 
sity and pressure are given by 



it 



Ao 



2n-4 



nuJo{2 — 2S ■ 



n 



2S + S^ + Q 



and 



(22) 



+ 2Sn-Q- ojan' 



45 - 52 - 3Q - n{6n -6 + 45) 



(29) 



These conditions give for p, p and 77 the following expres- 
sions: 



00 1" 



87r(2- 



nQ — n^S + nS — n^uo^S , (23) 



\--iSn - nS^ + 4S-S^-3Q + 3n 

247r(2 — n) 

—9n^ + 2Sn^ + 6n + iojon^ + Scoon^S — 6uJon^ 

S 



and 



167r 



1-5- 



+ jt (24) 



(25) 



We should note that the set of values (wo, n, Aq) will de- 
fine interesting cases, which we consider in the following. 



A. Einstein Kasner type universe 

Einstein theory is obtained for n — 0, ujq — > 
(j>Q — ^. In this case the solutions take the form 



t- 



16nG 



5^- 



r 



487rG' 



45 - 52 - 3g 



5 



and 



16nG 



1-5 



00 and 



(26) 



(27) 



(28) 



The corresponding particular cases are studied in de- 
tail in ref. We have shown above that the shear 
coefficient of viscosity is always negative. As we saw 
earlier, the DEC imply that -p + 0(, < p < p + 6^,. 
Then, from (^6|) , (|27| ) and ( p8| ) for p < p + 0^ we have 
that 1 — 5 < and from —p + 0^ < p we obtain that 
— 5^ — 25 + 3Q < 0. Thus, from these expressions it is 
easy to see that one always gets that rj < and p > 0. 



-3uJon^ , (30) 

respectively. 

From eqs. (|2^) and (|3^) we see that a barotropic state 
equation is satisfied [p — jp), and then 

_ 45 - 52 - 3Q - n(6n - 6 -I- 45) - Sluotl^ 

3(52 + 2Sn-Q- ujon^) ^ 

and the constraint for the constants is 



2nu;o(l - 5 - n/2) = 25 - 5^ - Q, 



(32) 



which is obtained from eq. (^2|). This solution must sat- 
isfy the second law of thermodynamics (77 > 0) and the 
DEC {p > and — 1 < 7 < 1), i.e. there are four con- 
ditions to be satisfied. It is easy to show that the three 
inequalities of the DEC are linearly dependent. Then we 
get from p > 



+ 2Sn-Q- ojon^ > 



(33) 



and from p < p 

452 - 45 + 67i2 - 671 + 10775 > 0. (34) 

Now, from the conditions 77 > and ( ^3| ) we obtain 

S{l + n)>Q + LUon^. (35) 

Thus, a physically consistent model must satisfy the con- 
ditions (H), (11), and (HI). 

It can be shown that the condition (|3^ ) is compatible 
with the inequality (|35|), but not with the condition ([34|). 
Therefore, in this particular model based on Brans-Dicke 
theory, there arc no physically reasonable solutions for 
a viscous Kasner type model. This is similar to what is 
found in Einstein's theory. 

C. Scalar-tensor Kasner type universe 

In the following we do not consider the bulk viscosity, 
since as is well known, the shear viscosity at early time is 
much greater than the bulk viscosity |L1| ; therefore, we 
could take ^ = 0. In this case we can write a barotropic 
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equation of state p = 7p, where, in agreement with the 
DEC, 7 lies in the interval — 1 < 7 < 1. 

It is easy to check that in this case there is a family of 
solutions which satisfy the second law of thermodynamics 
together with the DEC, simultaneously. For instance, we 
have found the following particular solutions. When the 
parameters u and n take the values luo = 0.1 and n = 
—0.1, it is found that the parameters pi — 0.1, p2 — 0.12 
and P3 = 0.8 represent a solution of the field equations. 
Here, all the relevant quantities, such that the effective 
pressure p, the shear viscosity r/, the energy density p, 
etc. all these quantities are physically acceptable. Thus, 
there are particular situations in scalar-tensor theories of 
gravity where are not found the problems appearing in 
Einstein theory. 

Note that if 77 = 0, i.e. n = 1 — S, one has the solution 

X (2^05* - t^o - ^05"^ + 25 - 52 - Q) , (36) 



from which wc get the barotropic equation of state, p — 
p, with — 1<1/S'<1- A dust dominated universe 
is obtained when S — > 00. If S* = 1 (n = 0), then the 
model becomes stiff matter dominated. In the latter case 
we obtain 



P = P 



(1-Q), 



(37) 



167rt2 

which coincides with equation (21) of ref. 

V. GENERATION OF ENTROPY 

In the following we shall determine the generation of 
entropy in a scalar-tensor theory. Here we include the 
shear coefficient of viscosity only, since ^ <C 77. 

It is well known that the production of entropy could 
be related to the anisotropy of the universe . In order 
to see this we introduce the entropy current four-vector 
as 



(38) 



where, as before, m'^ represents the four-velocity, the 
baryon number density, kg is the Boltzmann constant 
and a the nondimensional entropy per baryon. Since 



(^ = 0), it is obtained that 
where A = — ^ — 1. 



T]A, 



(39) 



(40) 



We note that from ( p5[ ) one has rjt^~^ = const. This 
fact permits us to write all the quantities if we know the 
value of 77 at a given moment of time. We use here the 
values given in ref [|j for some parameters at t = 1000 
s (the universe is characterized by ionized H and He in 
approximate equilibrium with radiation). At that time, 

the number densities of electrons and protons, Up — 

^ ' 1000 

10^^ CTTi^-^, the temperature T^^^^ ~ 4 x 10® K 



^^1000 

and 77i„„„ 
time t 



2.8 X 10^^ gem ^. Then we get at any 



,'1000\^" 



(41) 



Note that from (Eoh we have — = const. Thus 

I V ^ 

the generation of entropy at any time can be written as 



a{t) 



/lOOO 

n,{t) T{t) \~r 



n , T 
'^1000 1000 



3— n 



(42) 



In order to determine how the generation of entropy 
decreases as time goes on, we determine this generation 
of entropy for two different instants in the evolution of 
the universe. We consider the plasma era (radiation era), 
where the temperature ranges 10^" K >T > 4000 K pi| . 
For the recombination of hydrogen we have i,,^^ ~ 0. 126 x 



4 X 10 CTTT. the temperature 



1014, 

T,^^ ~ 4000 K. Then, for the ratio between values of 
generation of entropy at time of recombination t^^^ and 
the time t ~ 1000 we get for the model studied above 



(n = -0.1), that 



10 , which implies a drastic 



reduction in the rate of entropy production during the 
plasma era. 



VI. CONCLUSIONS 

We have studied the scalar-tensor theories for Kasner 
type metrics in the presence of a viscous fluid. We have 
explored the possibility of describing reasonable physi- 
cal models for a Kasner type universe when we keep the 
thermodynamics (r/ > and ^ > 0) together with the 
dominant energy condition in theories more general than 
Einstein's gravity. We concluded, for the specific model 
treated here, that in the Brans-Dicke theory, like in the 
Einstein theory, it is not possible to describe the growth 
of entropy. 

In the case of the scalar-tensor theory we have found, 
for our particular choice of the scalar field (p, the coupling 
function Lu{(j)) ~ ujq and the potential function A(0) = 
Ao(/)~^/", that it is possible to keep the thermodynamic 
conditions (»7 > and > 0) together with the DEC, 
as was described in Sect. IV. C. For these models we have 
found that, during the plasma era, i.e. from T ~ 10^° 
K to Tree ~ 4.000 K, a drastic reduction of the rate of 
entropy occurred. 
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